Emergence of U(l) symmetry in the 3D XY model with Zq anisotropy 
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We study the three-dimensional XY model with a Z q anisotropic term. At temperatures T < T c 
this dangerously irrelevant perturbation is relevant only above a length scale A, which diverges as a 
power of the correlation length; A ~ £ aq . Below A the order parameter is U(l) symmetric. We derive 
the full scaling function controlling the emergence of U(l) symmetry and use Monte Carlo results 
to extract the exponent a q for q = 4, . . . , 8. We find that a q sa a4,{q/A) 2 , with 04 only marginally 
larger than 1. We discuss these results in the context of U(l) symmetry at "deconfined" quantum 
critical points separating antiferromagnetic and valence-bond-solid states in quantum spin systems. 
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A salient feature of the recently proposed theory of 
"deconfined" quantum critical points, which separate 
Neel and valence-bond-solid (VBS) ground states of an- 
tiferromagnets on the square lattice, is the emergence 
of U(l) symmetry [l|. The VBS is either dimerized on 
columns or forms a square pattern with plaquettes of 
four strongly entangled spins [1, 0] ■ I n both cases there 
are four degenerate patterns and, thus, Z4 symmetry is 
broken. However, as the critical point is approached the 
theory predicts a length scale A, diverging faster than the 
correlation length, A ~ £°, a > 1, below which the dis- 
tinction between columnar and plaquette VBS states dis- 
appears. The nature of the VBS state is manifested only 
when coarse-graining on length-scales I > A, whereas for 
I < A the Z4 symmetry is unbroken and is replaced by an 
emergent U(l) symmetry characterizing the fluctuations 
between columnar and plaquette order. 

Quantum Monte Carlo simulations ^ of an S = 1/2 
Heisenberg model with four-spin couplings have recently 
provided concrete evidence for a continuous Neel-VBS 
transition, and also detected U(l) symmetry in the VBS 
order-parameter distribution P(D x ,D y ), where D x and 
D y are VBS order parameters for horizontal and vertical 
dimers. There is no trace of the expected Z4 anisotropy 
in the VBS phase — the distribution is ring shaped — 
although the finite-size scaling of the squared order pa- 
rameter shows that the system is long-range ordered. 
This can be interpreted as the largest studied lattice size 
L = 32 < A. A ring-shaped distribution was also found 
in simulations of an SU(N) generalization of the S — 1/2 
Heisenberg model [5[ — possibly a consequence of proxim- 
ity of this system to a deconfined quantum-critical point. 

In order to better understand the U(l) features of these 
VBS states, and to guide future studies of them, we here 
exploit a classical analogy. In the three-dimensional XY 
model including a ^-anisotropic term, 

H = -J^2cos(9 i -6 j )-h^2cQs(q6 i ), (1) 

(i,3) 1 

the anisotropy is dangerously irrelevant for q > 4 0, 0, H, 



B Ell; i-e., the universality class is that of the isotropic 
XY model but the perturbation is relevant for T < T c 
above a length-scale A. In the closely related g-state clock 
model, the anisotropy is dangerously irrelevant for q > 
5. While numerical studies [ll|, [l2|, EH have confirmed 
the irrelevance of the anisotropy at T c , the associated A 
has, to our knowledge, not been extracted numerically, 
except for an analysis of the 3-state antiferromagnetic 
Potts model, which corresponds to Zq [j| [3| • 

Here we report results of Monte Carlo simulations for 
4 < q < 8 on periodic-boundary lattices with N = L 3 
sites and L up to 32. In addition to Metroplis single-spin 
updates, we also use Wolff cluster updates [lj| to reduce 
critical slowing down. We sample the order-parameter 
distribution P(m x ,m y ), where 
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— cos(0i), m y =jjJ2 sin ^ 



(2) 
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The standard order parameter can be defined as 

f 1 f 1 ] 

(m) = / dm x / dm y P(m x ,m y ) (m x + m w ) 

dr I ddr 2 P(r,9). 
/o Jo 
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We will compare this with an order parameter (m q ) which 
is sensitive to the angular distribution; 
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= I dr I d9r z P(r,6)cos(q9). 



(4) 



While the finite-size scaling of (m) is governed by the 
correlation length £, (m 9 )should instead be controlled by 
the U(l) length scale A [9(, becoming large for a system 
of size L only when L > A. Fig. Q] shows magnetization 
histograms at h/ J = 1 for Z4 and Z§ systems with L = 4 
and 32. The angular distribution P{6) = J drrP(r, 9) 
is also shown. The average radius of the distribution is 
the magnetization (m), which decreases with increasing 
L. The anisotropy, on the other hand, increases with 
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FIG. 1: (Color online) P(m x ,m y ) at h/J = 1 for g = 4,8, 
L = 4,32. The temperature T/J = 2.17 for Z 4 and 1.15 for 
Zg; both less than T c /J ~ 2.20. The size of the histograms 
corresponds to m llS £ [—1,1]. Angular distributions P(0) 
with 6 € [0, 27r] are shown above each histogram. 

L. This is particularly striking for Zg, where the L = 4 
histogram shows essentially no angular dependence, even 
though T is very significantly below T c , whereas there 
are 8 prominent peaks for L = 32. Thus, in this case 
the U(l) length scale 4 < A < 32. For the Z 4 system T 
is much closer to T c but still some anisotropy is seen for 
L = 4; it becomes much more pronounced for L = 32. 

It is instructive to examine a spin configuration with 
m x « m y , i.e., 8 « n/4. Fig. [2] shows one layer of a Z 4 
system with L = 10 below T c . The spins align predomi- 
nantly along 8 — and 8 = it/ 2, with only a few spins in 
the other two directions. Clearly there is some cluster- 
ing of spins pointing in the same direction — the system 
consists of two interpenetrating clusters. Essentially, the 
configuration corresponds to a size-limited domain wall 
between 8 = and 8 = n/4 magnetized states. 

Hove and Sudb0 studied the g-state clock model and 
performed a course graining at criticality [HI]. They 
found that the structure in the angular distribution di- 
minished with the size of the block spins for q > 5, as 
would be expected if the anisotropy is irrelevant. Here 
we want to quantify the length scale A at which the 
anisotropy becomes relevant for T < T c . Consider first 
what would happen in a course graining procedure for 
a single spin configuration of an infinite system in the 
ordered state very close to T c . The individual spins will 
of course exhibit q preferred directions, as is seen clearly 
in Fig. [21 i.e., there would be q peaks in the probability 
distribution of angles Constructing block spins of I 3 
spins, we would expect the angular dependence to first 
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FIG. 2: (Color online) Spins in one layer of the Z± model 
with L = 10 at h/J = I, T/J = 1.9 < T c . Here m x w m y , 
corresponding to 8 ~ 7r/4 in P(r,9). Arrows are color-coded 
according to the closest Z4 angle; nn/2, n = 0, 1, 2, 3. 

become less pronounced because of the averaging over 
spins pointing in different directions (again, as is seen in 
Fig. [2]) . Sufficiently close to T c we would expect the dis- 
tribution to approach flatness. However, since we are in 
an ordered state, one of the q preferred angles eventually 
has to become predominant, and thus one peak in the 
histogram will start to grow. This happens at I » A. We 
cannot simulate the infinite system and instead we carry 
out an analogous procedure as a function of the lattice 
size L, sampling a large number of configurations. We 
calculate the order parameters (m) and (m 9 ), defined in 
Eqs. (|3I4|) , and analyze them using 

(m) - L-'fitL 1 ^), (5) 
(m q ) = L- a g(tL^). (6) 

Here ([5]) is the standard finite-size ansatz with a = / 3/v , 
and the XY exponents are (3 w 0.348 and v rj 0.672 [171 ]. 
Eq. ([6]) is an intuitive generalization of |(SJ), which was 
proposed and used also in Ref. |§], but we can actually 
also derive the scaling function g(X) exactly. 

Let us consider the scaling behavior of the order- 
parameter distribution P{fh). It depends upon the sys- 
tem size L and the size of scaling operators that perturb 
the critical theory. Specifically, we consider the temper- 
ature deviation t = T c — T and the presumed irrelevant 
q-fold anisotropy strength h. By conventional scaling ar- 
guments, we expect 

P(m; L, t, h) = L a/2 P(L a m, tL 1/u ,H = hL 3 - A "), (7) 

where A q > 3 is the scaling dimension of the irrelevant 
anisotropy. The prefactor above is determined from nor- 
malization of the probability distribution. In the scaling 
regime, |t| <C 1,L 1, so H is small. When the first 
two arguments are 0(1), P can be well-approximated by 
taking H = [with "corrections to scaling" of O(H), i.e. 
suppressed by L 3 ~ Aq for a large system]. At H = 0, 
the distribution is fully XY symmetric, and the integral 
in Eq. ([4]) vanishes. Thus, in this regime (m q ) is small, 
0(H), and should be considered as arising from correc- 
tions to scaling. This simply reflects the irrelevance of 
the anisotropy at the critical point. 
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Because the anisotropy is dangerously irrelevant, a 
larger contribution, however, emerges when tL x l v > 1, 
i.e. L ^> £ ~ t~ u . In this limit, the system can be re- 
garded as possessing long-range XY order, and the only 
significant fluctuations are the global fluctuations of the 
XY phase 9. This is biased by the anisotropy. The scale 
k of the total anisotropy (free) energy can be estimated 
by its typical magnitude within an XY correlation vol- 
ume, ft.£ 3_Ag , multiplied by the number of correlation 



volumes, (£/£) , i.e. K — hL < £ 



Note that although 



the energy per correlation volume is small (due to the ir- 
relevance of anisotropy at the critical point), the number 
of correlation volumes becomes very large and more than 
compensates for this smallness for L/£ sufficiently large. 

From this argument, we see that for L/£ 3> 1, the dis- 
tribution of angles is just determined from a Boltzmann 
factor for a single XY spin with the q-fold anisotropy 
energy ~ —ncosqO. Furthermore, for L/£ 3> 1, the 
magnitude \m\ ~ (to) is approximately non-fluctuating. 
Thus the distribution factors into the form P(\m\,9) = 
(m)~ 2 5{\m\ - (m))P(0), with 
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Here Z = J Q 71 dO e KCOS (i e ) [ s the single-spin partition func- 
tion. It is then straightforward to obtain from Eq. ([4]) 
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where /„ is the modified Bessel function of order n. Os- 
hokawa obtained a similar expression in a different way, 
but we disagree with his scaling variable. Comparing this 
with the scaling form in Eq. ©, we see that v q — A 9 z//3 
(a q = A g /3), k = hitl 1 /"*) 3 "* , and 



g(X) cx 



(10) 



Here h should be viewed as a non-universal scale factor. 
From the above discussion, one sees that this form is 
valid for L/£ ^> 1 but L/£ q arbitrary. For L/£ of O(l) or 
smaller, L/^ q <C 1 (implying n,X <C 1), and the scaling 
form for (m q ) becomes small and of order the expected 
correction to scaling in the critical regime. 

In Fig. [3] we show results for the two order parame- 
ters for systems with q — 4,5,6. We have studied several 
values of ft,/ J and here show results for a different value 
for each q. We have extracted T c using finite-size scaling 
of (m) with Eq. §5§ and the XY exponents. This works 
very well for all q, confirming the irrelevance of h. The 
magnetization for T < T c is seen to decrease marginally 
with increasing q in Fig.[3l The Z q order-parameter (m q ) 
changes more drastically, being strongly suppressed close 
to T c for large q. This is expected, as (m q ) should van- 
ish for all T in the XY limit q — > oo. For Z±, the (m q ) 
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FIG. 3: (Color online) The XY order parameter (m) (solid 
curves) and the Z q order parameter (m q ) (dashed curves) 
vs temperature for q = 4, 5, 6. The system sizes are L = 
8, 10, 12, 14, 16, 24, and 32. The curves become sharper (in- 
creasing slope) around T c (indicated by vertical lines). The 
ratios h/J used are indicated on the graphs. 



curves for different L cross each other, with the crossing 
points moving closer to T c as L increases. This is consis- 
tent with the above discussion of course-graining: In the 
ordered state close to T c , (m q ) should first, for small L, 
decrease with increasing L as the g-peaked structure in 
P{9) diminishes due to averaging over more spins. For 
larger L, (m q ) starts to grow with L as the length-scale 
A is exceeded. This behavior is more difficult to observe 
directly for q — 5, 6 because (m q ) is small and dominated 
by statistical noise close to T c where the curves cross. 

Fig. [1] shows finite-size scaling of the Z q order- 
parameter (m q ), using the hypothesis ^ and the XY 
value for a. Adjusting v q = a q v for each q we find sat- 
isfactory data collapse using 04 = 1.07(3), 05 = 1.6(1), 
a6 = 2.4(1), and, not shown in the figure, a§ = 4.2(3). 



These results are consistent with the form a q — 04(3/4) , 
in qualitative agreement with the e-expansion by Os- 
hikawa, which gave a q — > q 2 /30 for large q 0]. However, 
in the e-expansion there are significant deviations from 
the q 2 form in the range of q values considered here. Our 
ae is also smaller than the value ~ 3.5 obtained on the 
basis of the 3-state Potts antiferromagnet 

In Fig.[l]we also show the scaling function (fTO"]) . It does 
not match exactly the collapsed data, but the agreement 
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FIG. 4: (Color online) Scaling of the Z q magnetization. We 
use g = 0.52 for all q, and v q as indicated in the plots. The 
colors of the curves correspond to L as in Fig. [3] The dashed 
curves are the predicted scaling functions with h and prefac- 
tors adjusted to fit the data approximately. 



a VBS. Either a columnar or plaquette VBS should ob- 
tain in the infinite-size limit, but close to a deconfined 
quantum-critical point, for L < A, the system fluctuates 
among all mixtures of plaquette and columnar states. 
This corresponds to a ring-shaped VBS order-parameter 
histogram. In numerical studies of quantum antiferro- 
magnets 0, [H| no 4-peak structure was observed in the 
angular distribution, and hence it is not clear what type 
of VBS finally will emerge (although a method using open 
boundaries favors a columnar state in fH). It seems un- 
likely that the U(l) symmetry should persist as L — > oo. 
In the classical Z$ model we never observe a perfectly 
U(l)-symmetric histograms far inside the ordered phase, 
in contrast to Refs. [4j,|5[. On the other hand, a q is larger 
for q > 4, and in Fig. [T] we have shown a prominently 
U(l)-symmetric histogram for the Zg, model deep inside 
the ordered phase. Thus, the exponent a may be larger 
for the Z4 quantum VBS than 04 ~ 1 obtained here for 
the classical Z4 model. There is of course no reason to 
expect them to be the same, as the universality class of 
deconfined quantum-criticality is not that of the classical 
Z4 model [lj, |4|. Future numerical studies of VBS states 
and deconfined quantum-criticality can hopefully reach 
sufficiently large lattices to extract the U(l) exponent 
using the scaling method employed here. 
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improves as q increases. As we have discussed above, 
the scaling function represents the dominant behavior 
for T < T c , but exactly at T c this contribution vanishes 
and the critical-point scaling form (m q (T c )) ~ L 3 ~ Aq be- 
comes dominant. For q — 4, A4 — 3 is small; our estimate 
is A4 — 3 = 0.21(9), in good agreement with previous es- 
timates of the scaling dimension 0, [l(| • Thus it is clear 
that very large systems would be required for this contri- 
bution to become invisible on the scale used in our graph. 
It is also clear that for t < 0, close to t = 0, there will 
be a similarly significant correction to the asymptotically 
dominant scaling function. As q increases, we have seen 
that A q = 3a q increases rapidly, and we thus expect sig- 
nificantly smaller correction to scaling. For q = 6 the 
agreement is already seen to be quite good, considering 
that our lattices are not very large. 

To conclude, we relate our results to the quantum VBS 
states discussed in the introduction. Returning to Fig. [21 
associating rs arrows with two adjacent horizontal 
dimers on even- numbered columns and 0i m n/2 with 
vertical adjacent dimers on even rows, (0) = 0,7r/2 cor- 
respond to columnar VBS states. A plaquette is a super- 
position of horizontal and vertical dimer pairs, whence a 
plaquette VBS corresponds to (6) — 7r/4 3]. Rotating 
the arrows by 90° corresponds to translating or rotating 
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